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^ ' Abstract 

(N 
(N 

(^ . We give bounds for the global attractor of the delay differential equation x{t) = —^x{t) + 

CO ' f {x{t—T)) , where / is unimodal and has negative Schwarzian derivative. If / and jj. satisfy certain 

r>~^ . condition, then, regardless of the delay, all solutions enter the domain where / is monotone 

(^ ' decreasing and the powerful results for delayed monotone feedback can be applied to describe 

OO , the asymptotic behaviour of solutions. In this situation we determine the sharpest interval that 

contains the global attractor for any delay. In the absence of that condition, improving earlier 
results, we show that if the delay is sufficiently small, then all solution enter the domain where /' 
is negative. Our theorems then are illustrated by numerical examples using Nicholson's blowflies 
equation and the Mackey- Glass equation. 
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1 Introduction 

This note is motivated by a recent paper by G. Rost and J. Wu [IT] about the so-called delayed 
recruitment model defined by the delay differential equation 

x'(i) = -/ix(t) + /(x(t-r)), (1.1) 

where /U > 0, r > 0, and / : [0, oo) -^ [0, cx)) is a continuous function. 

In particular, they consider the case when / is a unimodal function, which is the situation for the 
famous Nicholson's blowflies equation and the Mackey-Glass model. 



In that reference, the authors have proved several results on the global dynamics of Eq. (jl.ip . and 
they also formulated some open problems. It is our purpose to prove new results in the direction 
initiated in [11], and also to answer some of the open questions. We show the applicability of our 
results for different cases of the Nicholson's blowflies equation 

x'{t) = -nx{t) +px{t- r)e-^^(*-^) , (1.2) 

where fi, p, 7 are positive parameters (see, e.g., [7] for a biological interpretation); and the Mackey- 
Glass equation [8] 

x'{t) = -f^xit) + ^f^;~^\„ . (1.3) 

1 + X[t — Tj" 

Following [H], we assume that / is unimodal. More precisely, the following hypothesis will be 
required: 

(U) f{x) > for all x > 0, /(O) = 0, and there is a unique xq > such that f'{x) > if 
< X < xo, /'(xo) = 0, and /'(x) < if x > xq- Moreover, /"(x) < if < x < xq. 

Much is known about the global picture of the dynamics of Eq. (jl.ip . when / is a monotone 
function. However, unimodal feedback may lead to very complicated and still not completely 
understood dynamics. See [H [5l [6] and references thereof. 

We shall use the function g(x) = fj.~^f{x). Notice that the equilibria of (jl.ip are the fixed points 
of g, and that, under condition (U), function g has at most two fixed points x = and x = K > 0. 

As in [11], we consider nonnegative solutions of (jl.ip . We recall that for each nonnegative and 
nonzero function (/> G C = C([— r, 0],M), there exists a unique solution x'^{t) of (jl.ip such that 
x<^ = ,/. on [-r,0]. Moreover, x<^(t) > 0, Vi > 0. (See, e.g., \T, Corollary 12].) 

It is well known (see, e. g., [2l[TT]) that all solutions of (jl.ip converge to if 5f'(0) < 1, whereas the 
positive equilibrium K is globally attracting for Eq. (jl.ip if g'{K) > (equivalently, if i^ < xq). 

Thus, we shall assume that g'{0) > 1 and K > xq. Ivanov and Sharkovsky [Sj Theorem 2.3] proved 
that an invariant and attracting interval [a, (3] for g is also invariant and attracting for (jl.ip for all 
values of the delay r, that is, 

a < liminf x(t) < limsupx(i) < /?, 



t— >oo 



i— >oo 



for any nonzero solution x of (jl.ip . Similar results were proven using a slightly different approach 
in [21 E]. 

It is clear that we can choose /? = g{xo), a = g{l3) = g'^{xo) to get an attracting invariant interval 
[a, (3] for the map g (here, and in the following, g^ denotes the composition gog). Thus, this interval 
contains the global attractor associated to Eq (jl.ip for all values of r. Using this fact, Rost and 
Wu obtain sufficient conditions to ensure that every solution of (jl.ip enters the domain where /' is 
negative. In this case the asymptotic behaviour of the solutions is governed by monotone delayed 



feedback and the comprehensive theory of monotone dynamics is apphcable, as it was demonstrated 
in [11]. In particular, since a Poincare-Bendixson type theorem is available for (jl.ip when /' is 
negative, this kind of conditions exclude the possibility of solutions with complicated asymptotic 
behaviour (the w-limit set can only be the positive equilibrium K oi a periodic orbit). We include 
here the main results of [IT] in this direction 

Theorem 1. Every solution of lil-l]) enters the domain where f is negative if any of the following 
conditions holds: 



(L) a = g{[3) = g^{xo) > xq. 
n(xo) - xq 



where H : {0,K) — > [K,oo) is the inverse of the restriction of g to the interval [K,oo). 

Notice that the first condition in Theorem [1] is independent of the delay, while condition (L,-) shows 
that even if / is unimodal, the solutions of (jl.ip have the same asymptotic behaviour as in the case 
of monotone decreasing feedback for all sufficiently small delay r. 

An open problem suggested in [11] is the following: under condition (L), find the sharpest invariant 
and attracting interval containing the global attractor of (jl.ip for all r. (Numerical experiments 
performed in [llj show that J = [a,/3] seems to be a very sharp bound). To avoid confusion we 
remark that the global attractor ^ is a subset of the function space C([— r, 0],M), and saying that 
an interval [a, b] contains the global attractor we mean that for each i;^ G ^, we have a < 0(s) < b 
for any s € [— r, 0]. 

Our main results in this note are the following: 



1. We completely solve this problem in the case of Nicholson's blowflies and Mackey-Glass equa- 
tions, obtaining an interesting dichotomy result for (II. 2p and (11. 3p when condition (L) holds 
(Theorem [6]). 

2. We give a weaker delay-dependent condition different from (L,-) under which the statement 
of Theorem [1] remains valid. In other words, we can determine a r* that is larger than r* in 
Theorem 1, such that all solutions enter the domain where /' is negative, if r < r*. Moreover, 
we provide some examples showing that the new condition significantly improves (L,-) in 
certain situations. 



2 Main Results 



In this section, we assume that condition (U) holds, ^'(0) > 1 and K > xq, where g = fj, ^ f . 
Denote 

a = ini{A>Q : g^{A) = A} ; ^ = sviv{B > Q : g^{B) = B}. (2.1) 



Since ^'(0) > 1, a and (3 are well defined real numbers. Assuming that (L) is satisfied, it is clear 
that a > xq and hence g is decreasing on J := [a,/3]. Moreover, ^(q) = /?, and g{(3) = a. As a 
consequence, J is invariant for the map g. 

Lemma 2. Assume that (L) holds. Then, J = [a, (3] is an attracting invariant interval for the 
map g. 

Proof. We have already proved that J is invariant. Next, we prove that it is attracting. 

Since h = g'^ is monotone increasing in J, and a and f3 are respectively the minimal and the 
maximal fixed points of /i in [a,/?], it follows that 

lim /i"(x) = a, Vx G [a, a] ; lim /i"(x) = ^5, Vrc G [13, /S]. 

n—>oo n— ►oo 

Since g{[a, a]) = [(3, (3] and g{[(3, (3]) = [a, a], the result follows from the fact that [a, /3] is attracting 
for g. n 

An application of the above mentioned Theorem 2.3 in [3] gives the following result: 

Corollary 3. Under condition (L), the interval J = [a, /3] is an attracting and invariant interval 
for M.l\) for all values of the delay r. 

Remark 1. If the equilibrium K is globally attracting for g, then a = f3 = K , and hence K is also 
a global attractor for lil.l]) for all values of the delay. This result is |3, Theorem 2.2]. 

Remark 2. // (L) does not hold, the interval J does not need to be globally attracting for lil.l\l . 
For example, for the Nicholson's blowflies equation considered in flT^ 

x'{t) = -0.052;(t) + x{t - T)e-<^-^\ (2.2) 

the interval J is given by 

J = [a,/3] ^ [0.4261,5.5653]. 

It is easy to check that this interval does not attract every orbit associated to the map g, since there 
is a period four orbit given by 

0.24286 -^ 3.80991 -^ 1.6878 -^ 6.24235 -^ 0.24286. 

Numerical experiments from [TTf suggest that J also does not attract every orbit of i2.^) for large 
values of T. 

We recall here that the nonlinearity / in some important examples of Eq. (II. ip (including the 
Mackey-Glass and Nicholson's blowfiies models) fulfills the following additional assumption: 



(S) / is three times differentiable, and {Sf){x) < whenever f'{x) / 0, where Sf denotes the 
Schwarzian derivative of /, defined by 



fix) 2V/'W 



The following proposition is a consequence of Singer's results [T2] . 

Proposition 4. Assume that g : [a,b] -^ [a,b] satisfies (S) and g'{x) < for all x G [«,&]■ Let K 
he the unique fixed point of g in [a, h] . Then, 

• If \g'{K)\ < I, then Imin^oo g"^ {x) = K , M x (^ [a,b]. 

• If \g'{K)\ > 1, then there exists a globally attracting 2— cycle {p,q}. More precisely, g{p) = q, 

9{q) =P, Pt^Q, and 

( p if X < K 
lim g^'^ix) = { 
"-"°° [ q if x> K. 

For a detailed proof of the first statement in a more general situation, see, e.g., [71, Proposition 3.3]. 
The second statement can be easily proved using the same arguments. 

Remark 3. As noticed above, condition (S) holds for the Mackey-Glass and Nicholson's blowflies 
models, among others. See, for example, m [^ IJOJ/ . Thus, in these models, the interval J given in 
LemmalE reduces to {K} if\g'{K)\ < 1, whereas J = [p,q] if\g'{K)\ > 1. 

The next result shows that J is actually the sharpest invariant and attracting interval containing 
the global attractor of (jl.ip for all r if (L) and (S) hold. We emphasize that in the limit case 
5^(xo) = xq in (L), we have a = a = xq, (3 = (3 = g{xQ). Thus, as suggested in [11], the intervals J 
and J coincide in this special situation. 

Proposition 5. Assume that (L) is fulfilled, and (S) holds in the interval [a,0\. If\g'{K)\ > 1 
then, for any ^ > 0, there exists a sufficiently large rg such that the interval [a + C, /? — C] is not 
attracting for lil.l\) if t > r^ . 

Proof. We make use of Theorems 2.1 and 2.2 in flO] (see also [9]). Define e = (^r)~^. A direct 
application of [ini Theorem 2.1] provides an eq > such that for each e < Eq, Eq. (jl.ip possesses 
a slowly oscillating periodic solution p{t) and there exist constants C,'^ > satisfying p{t) G J, for 
all t, p{0) = p{C) = p{oj) = K, p{t) > K in (0, C), p{t) < K m. (C, w), p{t) = p(t + to) for ah t. 

Next, [101 Theorem 2.2] states that, given any 5 > 0, there exist ei > 0, A: > such that \p{t)—a\ < 6 
in [£k,C — sk] for all e < £i. Choosing 5 = ^/2, it is clear that the periodic solution p{t) is not 
attracted by the interval [a + ^, /3 — ^] if r > r^ = {fiei)~^. D 

Remark 4. What is important in the proof of Proposition is the fact that g has a unique glob- 
ally attracting 2-periodic solution. Condition (S) implies this fact, although it is not a necessary 
condition. For example, function 

qix) = 30 ^ 

satisfies (L) and (U), \g'{K)\ > 1, and it has a unique globally attracting 2-cycle defined by a = 
0.728449, /3 = 2.2822. The conclusion of Proposition \^ holds although {Sg){x) > for x > 2.02. 
Notice that [a,/5] = [0.515162,3.62133]. 



On the other hand, if g does not have a unique globally attracting 2-periodic solution, under con- 
dition (L), interval J is still the smallest globally attracting interval for the difference equation 
Xn+i = gixn), n = 0,1, . . . We conjecture that the conclusion of Proposition^^ remains valid with- 
out assumption (S). 

As a consequence of Corollary[3l Reniark[T]and Propositions[l]and[5l we get the following dichotomy 
for Equation (|l.ip under conditions (L) and (S): 

Theorem 6. Assume that (L) is fulfilled and (S) holds in the interval [a, /?]. Then exactly one of 
the following holds: 

(1) \g'{K)\ < 1 and the global attractor of l\l.l\) for all values of the delay r is {K}. 

(2) \g'{K)\ > 1 and the sharpest invariant and attracting interval containing the global attractor 
of UnP for all values of the delay t is [a, (3], where {a,/3} is the unique 2-cycle of g in [a, (3]. 

According to Theorem [H when condition (L) does not hold, it is still possible to find a delay- 
dependent condition (L^) under which every solution of (11. ip enters the domain where /' is negative 
and the theory of monotone delayed feedback can be applied to describe the asymptotic behaviour of 
our equation with unimodal feedback [11]. Next we give a different condition in the same direction. 
The proof is based on the following lemma proved in [2] (see also [TJ Lemma 5.1]). 

Lemma 7. Assume that (U) holds, g'{0) > 1, and K > xq. Then, for every nonnegative and 
nonzero solution x{t) of ( fj. Jl) there exist finite positive limits 

M = limsupa;(t) , m = liminf x(t). 



t^oo 



i— >oo 



Moreover, [m,M] C gi{[m,M]), where 

gi{x) := (1 - e-'^^) g{x) + e-^^^K. (2.3) 

Theorem 8. Assume that the following condition holds: 

(LV) glixo) > Xq, where gi{x) = (1 — e~^^) g{x) + e~^^'' K. 

Then, every solution of \1.1\) enters the domain where f is negative. 

Proof. Let x{t) be a nonnegative nonzero solution of (jl.ip . m = liminf x(t) and M = limsupx(t). 



i— >oo 



t— +00 



Notice that gi{x) lies between g{x) and K. In particular, 51 (xq) > xq. On the other hand, since 
/ satisfies (U), it is clear that gi also meets the same condition, except that 51 (0) > 0. Moreover, 
gi has the only fixed point K. Since, by Lemma[71 [m,M] C gi{[m,M]), we obtain m < K < M, 
thus m < gi{m), M > gi{M) and the following inequalities must hold: 

vcAn{gi{M),gl{xo)}<m, M < gi{xo). 




g'(xo) 



Xo K n(Xo) g(xj 



Figure 1: The profile of the function 5 in a situation where condition (L) does not hold. 

Since gi is decreasing in [M,gi{xQ)], M < gi{xo) implies gl{xQ) < gi{M) and then 

glixo) = m.m.{gi{M),gl{xo)} < m. 

Finally, (LV) implies that m > gf^xo) > xq, and therefore x{t) > xq for all sufficiently large t. The 
proof is complete. D 

As a byproduct of the proof of Theorem [8l we get the following result, which is of independent 
interest to obtain sharper bounds for the global attractor of Eq. (jl.ip . 

Corollary 9. Let g\ he the function defined in i2.3\) . If (U) holds, g'{0) > 1, and K > xq, then 
the interval [gi{xo),gi{xo)] contains the global attractor of lll.l\) . 

Notice that (LV) always holds if (L) is satisfied, since g'f(xo) > a = g^(xo). Next we show that 
when (L) does not hold, (LV) is sharper than (L,-), and therefore Theorem [8] improves Theorem 

3.8of [llj. 

Proposition 10. Assume that (L) does not hold. Then (Lr) implies (L'^)- In other words, (L'^J 
gives a better estimate than (Lr) for the possible delays that still guarantee that every solution 
enters the domain where f is negative. 

Proof. If (L) fails, then fl'^(xo) < xq. Thus we have g'^{xo) < xq < K < n(xo) < g{xo) (as depicted 
in Figure 1). Using the notation 6 = fir > 0, (L,-) is equivalent with 

Xq + 9{g{xo) - g^{xo)) < n(xo) 

and {L'^) can be written as 

(1 - e-')g[{l - e-')g{xo) + e'' K] + e-' K > xq. 



First we suppose that (1-e ^)g{xo) + e ^K < n(xo). Notice that ^(xo) > (1-e ^)g{xo) + e ^K > 
K, and g is decreasing in in [K,g{xo)], so 

g[il - e-')gixo) + e-'K] > g(n(xo)) = xq, 



and by 

(1 - e-^)g[{l - e-'^)g{xo) + e-^ K] + e-^ K > (1 - e-^)xQ + e'^ K > xq, 

{L'^) holds. So far we have not used (Lt), (1 — e~^)g{xo) + e~^K < n(xo) always implies {L'^)- 

Now we consider the remaining case (1 — e~^)g(xQ) + e~^K > n(xo). Since (L,-) holds, we have 
(1 - e~^)g{xo) + e~^K > xq + 0{g{xQ) - g'^{xo)). On the other hand, 

g[il-e-')g{xo) + e-'K]>g\xo). 

Therefore, 

(1 - e-')g[{l - e-')g{xo) + e'' K] + e'' K > (1 - e-')g\xo) + e-' K = 
= (1 - e-')g{xo) + e-'K + (1 - e-'){g\xo) - g{xo)) > 
>xo + e{g{xQ) - g\xo)) + (1 - e-^){g\xo) - g{xQ)) = 
= xo + (gixo) - g^ixo)){e + e-^ - 1) > xo, 

where in the last step we used that g^xo) — ^^(xo) > (this follows because (L) does not hold), 
and that the function h{6) = 9 + e^ — 1 is nonnegative for 6 > 0. D 

Remark 5. Notice that (L,-) always fails if > 1 (or equivalently r > l/fJ-), but we have not used 
this fact in the proof. 



3 Examples 

In this section, we use the Nicholson's blowflies equation and the Mackey-Glass equation with 
different parameters in order to illustrate our results in Section 2. 

After a change of variables, one can always write (II. 2p in the form 

x'(t) = -fix{t) + x{t - T)e~^(*-"). (3.1) 



It is well known that there is a unique equilibrium x = if ;U > 1, and it attracts all nonnegative 
solutions. If /x < 1, there is a positive equilibrium K = — ln(/i), and x = becomes unstable. 
Moreover, K is globally attracting for all values of the delay if ^ G (e~^, 1) ~ (0.13533, 1). Next, 
using Theorem 9.3 in flO], it follows that condition (L) holds for fj, S (i/, e~^), where v ~ 0.10472 
is defined by the relation v = e~^v~^, v being the unique solution greater than 2 of equation 
v^ = e^-i. 

For /i = 0.13, condition (L) holds, and the invariant and attracting interval for (j3.ip given by 
Theorem 3.5 in [llj is 

[a,/3] = [g\x^),g{xo)] = [1.2848,2.8298]. 



This interval may be strengthened, since the smallest invariant and attracting interval for (|3.1 
independent of the delay r is given by the unique 2-cycle of g, that is, 

[a,/3] ^ [1.54796,2.53248]. 



According to Proposition [5l there exists a slowly oscillating periodic solution of ()3.ip whose mini- 
mum and maximum values get closer and closer to a and f3 as t tends to infinity. See Figure 2. A, 
where two distinct solutions are presented and the horizontal lines indicate a and /?. 



A) Nicholson (t=100, |i=0.13) 




B) Mackey-Glass (t=4, |i=1.79) 




500 1000 1500 2000 



10 20 30 40 50 60 70 



Figure 2: Nicholson's blowflies equation and Mackey-Glass equation when the condition (L) holds 

For /x = 1/16 = 0.0625, condition (L) does not hold (indeed, g'^{xo) = 0.2616 < xq = 1). Since 
n(xo) = 4.21007, we get 

* _ n(xo) - xo 



fJ-igixo) -g^{xo)) 



0.570734. 



On the other hand, one can check that condition (LV) in Theorem [8] holds for t < t^, = 1.46534, 
showing how (LV) gives an estimate significantly sharper than (L^). 

It is interesting to notice that from Theorem 2.1 in [7] it follows that the positive equilibrium K 
is globally attracting for (j3.ip if r < 1.1935. Hence, the information provided by Theorem[T]is not 
very useful here, whereas Theorem [8] can be applied for r between 1.1935 and 1.46534. 

Our next example is the Mackey-Glass equation 



A) Mackey-Glass (t=0.195, h=1) 



B) Mackey-Glass (t=3, h=1) 





10 20 30 40 50 



Figure 3: Illustration of Theorem [8) Mackey-Glass equation with different delays 



2x{t - r) 
1 + xit- rf 



x'{t) = -/xx(t) + , _7, l!,^ . (3.2) 



The function f{x) = 2x(l + x^'^)^^ satisfies the unimodal condition (U) with xq = 0.863. For ^ > 2, 
there is a unique equihbrium x = which attracts all nonnegative solutions. If /i < 2, there is also 
a positive equilibrium K = (—1 + 2/2"^)^"^^. Next, K < xq for /i G [1.9, 2) (and hence K is globally 
attracting), while i^ > xq for ^ < 1.9. 

Denote, as usual, g = fi^^ f. One can check that (L) holds if /i > 1.774 so in the interval (1.774, 1.9) 
the dichotomy stated in Theorem [6] applies. Since |g''(ir)| < 1 if and only if /i > 1.8, the equilibrium 
K attracts all positive solutions for fi G [1.8, 1.9), whereas for each fi G (1.774, 1.8) we can determine 
the sharpest delay- independent interval [a, /3] containing the global attractor of (j3.2p by finding the 
unique 2-cycle of g in the interval [a,/3] = [g'^{xQ),g{xo)]. For example, setting /i = 1.79, we 
have K = 0.898, [a,f3] = [0.872,0.916] and the interval [a,/3] = [0.876,0.914] contains the global 
attractor of (j3.2p for all values of the delay. See Figure 2 B, where the horizontal lines represent a 
and p. 

When (L) does not hold, we still can use Theorem[8j For example, for /x = 1 the positive equilibrium 
K = 1 loses its asymptotic stability for r > 0.188. One can check that condition (LV) in Theorem[8] 
holds for r < n = 0.195, so for r G (0, n) every solution of (|3.2p enters the domain where /' is 
negative, while (L,-) is satisfied only for t < t* = 0.092. In this situation apparently chaotic 
behavior can be observed for large enough delays (see Figure 3 B with r = 3), but Theorem [8] 
and the results of [IT] guarantee that complicated behavior is not possible if r < r* = 0.195, 
as illustrated in Figure 3 A. Moreover, we get a good bound for the global attractor of (j3.2p 
from Corollary m Indeed, in this case the interval [a,f3] = [0.00016,1.639] is improved up to 
[gi{xo),gi{xo)] = [0.864,1.113]. See Figure 3 A, where the horizontal lines indicate qKxq) and 
9i{xo). 
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